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We develop a method for investigating the relationship between the shape of a 1-particle distribution 
^3 ' and non-linear electrostatic oscillations in a collisionless plasma, incorporating transverse thermal motion. 

C^ I A general expression is found for the maximum sustainable electric field, and is evaluated for a particular 

highly anisotropic distribution. 

Introduction 

C/3 , High-power lasers and plasmas may be used to accelerate electrons by electric fields that are orders of magnitude 
^ ' greater than those achievable using conventional methods [1]. An intense laser pulse is used to drive a wave 
Qh! in an underdense plasma and, for sufficiently large fields, non-linearities lead to collapse of the wave structure 
C/3 ' ("wave-breaking") due to sufficiently large numbers of electrons becoming trapped in the wave. 
O I Hydrodynamic investigations of wave-breaking were first undertaken for cold plasmas [2, 3] and thermal 

C/3 ' effects were later included in non-relativistic [4] and relativistic contexts [5-7] (see [8] for a discussion of the 
^~»| numerous approaches). However, it is clear that the value of the electric field at which the wave breaks (the 
electric field's "wave-breaking limit") is highly sensitive to the details of the hydrodynamic model. 

Plasmas dominated by collisions are described by a pressure tensor that does not deviate far from isotropy, 
whereas an intense and ultrashort laser pulse propagating through an underdense plasma will drive the plasma 
►^ ! anisotropically over typical acceleration timescales. Thus, it is important to accommodate 3-dimensionality and 
^—v ' allow for anisotropy when investigating wave-breaking limits. The sensitivity of the wave-breaking limit to the 
^^ I details of the plasma model suggests that it could depend on the anisotropy of the pressure tensor. 
t^^ ' One method for investigating the wave-breaking limit of a collisionless anisotropic plasma is to employ 

CO I the warm plasma closure of velocity moments of the 1-particle distribution / satisfying the Vlasov-Maxwell 
[■*-. ' equations [7] . Successive order moments of the Vlasov equation induce an infinite hierarchy of field equations 
^D I for the velocity moments of / and at each finite order the number of unknowns is greater than the number of 
field equations. The warm plasma closure scheme sets the number of unknowns equal to the number of field 
equations by assuming that the terms containing the third order centred moment are negligible relative to those 
including second, first and zeroth order centred moments. 

Our aim is to uncover the relationship between wave-breaking and the shape of /. In general, the detailed 

; ^ ' structure of / cannot be reconstructed from a few low-order moments so we adopt a different approach based on 

C^ I a particular class of piecewise constant 1-particle distributions. Our choice of distribution, although somewhat 

artificial, reduces the Vlasov equation to that of a boundary in the unit hyperboloid bundle over spacetime. 

Combining the equation for the boundary with the Maxwell equations yields an integral for the wave-breaking 

limit in terms of the shape of the boundary. 

Our approach may be considered as a multi-dimensional generalization of the 1-dimensional relativistic 
"waterbag" model employed in [5]. 

1 Vlasov-Maxwell equations 

The brief summary of the Vlasov-Maxwell equations given below establishes our conventions. Further discussion 
of relativistic kinetic theory may be found in, for example, [9,10] . We employ the Einstein summation convention 
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throughout and units are used in which the speed of hght c ~ 1 and the permittivity of the vacuum Eq = 1. 
Lowercase Latin indices a,b,c run over 0,1,2,3. 

Preliminary considerations 

Let (a;°) be an inertial coordinate system on Minkowski spacetime {Ai,g) where a;° is the proper time of 
observers at fixed Cartesian coordinates {x^ ,x'^,x^) in the laboratory. The metric tensor g has the form 

g = 7]abdx'' (g)dx'' (1) 

where 

-lifa = & = 
T^ab^ { liia = b^O (2) 

Oifa^b 

Let (x'^jx'^) be an induced coordinate system on the total space TA4 of the tangent bundle {TA4,Il,A4) 
and in the following, where convenient, we will write x instead of x" and x instead of x''. 

We are interested in the evolution of a thermal plasma over timescales during which the motion of the ions is 
negligible in comparison with the motion of the electrons. We assume that the ions are at rest and distributed 
homogeneously in the laboratory frame. Their worldlines are trajectories of the vector field Mon = n-mad/dx^ on 
M. where njon is the constant ion number density measured in the laboratory frame. The electrons are described 
statistically by a 1-particle distribution f{x,x) which induces a number 4-current vector field N = N°'d/dx°' 

7^(2;)= / x''f{x,x)^^^^=dx^dx^dx^ (3) 

on Ai, where \x\'^ = (.-i^)^ + (i^)^ + {x^Y . One may write the Maxwell equations on M as 
dFbc , dFab . dF,, 



ca 



dx'^ dx'^ dx^ 

Qpba 



-0, (4) 



dx' 



- qN'^ - qN,l^, (5) 



where Fab are the components of the electromagnetic field tensor, F"** ~ ri°''^Tf"^Fcd, Q is the charge on the 
electron (g < 0) and (77°'') is the matrix inverse of [rjab)- The scalar field / satisfies the Vlasov equation, which 
may be written 

\dx'' m "■ dx^ J ^ ' 

on TM where F^"^ is the vertical lift of F\ = if^F^a from M to TM, 

F'Y{x,x) = F\{x). (7) 

Exterior formulation 

In this section we recast the above using the tools of exterior differential calculus as it affords a succinct and 
powerful language for subsequent analysis. We make extensive use of Cartan's exterior derivative d, the exterior 
product A and the Hodge map * on differential forms (see, for example, [11, 12]). 
The spacetime volume 4-form •I is 

•1 = dx^ A dx^ A dx^ A dx^ (8) 

and the Maxwell equations ([H [5]) can be written 

dF = 0, d*F ^ -qirN + qirNZn (9) 

where F = ^Fab dx"^ A dx'' is the electromagnetic 2-form, and the 1-forms N, Nion are the metric duals of the 
vector fields N, A^ion respectively. (The metric dual y of a vector field Y satisfies Y{Z) ~ g{Y, Z) for all vector 
fields Z.) 



Introduce the vector fields L, X, 

on TM and the 6-form uj, 

u^LLixi*!"^ A#l) (12) 

on TA4 where iy is the interior operator on forms with respect to vector Y, the 4- form •l^ 

•1^ = dx° A dx^ A dx^ A dx^ (13) 

is the vertical lift of the spacetimc volume 4- form •! from Ai to TA4 and the 4- form #1 

#1 = dx" A dx^ A dx^ A di^ (14) 

onTA^. 

The total space £ of the sub-bundle (£,11, A^) of {TAi,Il,Ai) is the set of timelikc, future-directed, unit 
normalized tangent vectors on A4 , 

£ = {{x, x) e TM I (^ = and i° > 0} (15) 

where 

(^ = r/afci"i;'' + l. (16) 

The integral ([3]) can be written 

N'^{x)= f xyix#l (17) 

where £x = n^^(x) is the fibre of {£, 11, Al) over x G Al, and it can be shown that the Vlasov equation ^ can 
be written 

difio) ~ (18) 

where ~ denotes equality under restriction to £ by pull-back. Thus, it follows 

d{fLo) = (19) 

where ,S is a 6-dimensional region in £ and using the generalized Stokes theorem on forms (see, for example, [12]) 
we obtain 

fio^O (20) 

dB 

where dB is the boundary of B. 

Piecewise constant distributions 

We consider distributions for which / = a is a positive constant inside a 6-dimensional region hi d £ and / = 
outside. In particular, we consider U to be the union over each point a; G A4 of a domain Wx whose boundary 
dWx in £ is topologically equivalent to the 2-sphere. Such distributions are sometimes called "waterbags" in 
the literature. 

Choosing B in ([20]) to be a small 6-dimensional "pill-box" that intersects dWx and taking the appropriate 
limit as the volume of B tends to zero, we recover a jump condition on fuj that leads to 

d\Auj~Q atA = (21) 



where A = is the union over x of the boundaries dWx- 
If A = is the image of the embedding map S, 

S : TW X 52 ^ S 

(x,0 ^ ix,x = V^{x)), (22) 

where ^ = (^\^^) is a point in 5*^, tlien it follows from ([101 El [12|) that ([2l|) is equivalent to 

NvM-—i'V,F)An.=0. (23) 

' m ^ 

Here, V^ and fl^ are families of vector fields and 2-forms on Ai respectively, defined by 

y^ = yti^ = (^*^')i^^ (24) 

where dxa = rjabdx^ . Note that since the image of S lies in 8, it follows that, for each ^ e S"^, V^ is timelike, 
unit normalized and future-directed: 

.g(Fj,y«) = -l, 5(F^,A)<o. (26) 



We adopt ([23]) as the equation of motion for dWx ■ 

It may be shown that a particular class of solutions to ([^5)1 satisfies 



777 — 
F=-dV^ (27) 

and using Q we obtain the field equation 

„2 ^ 

di,dVe = - — {-*:N ^ *iVio„) (28) 

ra 

on A^ with the condition that dVj is independent of ^. For simplicity, wc have neglected the direct contribution 
of the laser pulse to the total electromagnetic field in ([77)1 . 

2 Electrostatic oscillations 

Before analysing (|281 I26p further it is useful to briefly discuss their analogue on 2-dinicnsional spacetimc for 
facilitating comparison with the approach adopted in [5]. 

Electrostatic oscillations in 1 spatial dimension 



Although formulated on 4-dimensional spacetime, equations pSl I26|) have a similar structure for any number 
of dimensions. In particular, we now consider 2-dimensional Minkowski spacetimc (A^,g) 

g = -dt ®dt + dz® dz, (29) 

•1 = dtAdz (30) 

where (t, z j^ is a Cartesian coordinate system in the laboratory inertial frame. An induced coordinate system 
on TA4 is (t, z, t, z) and note that in this sub-section of the article the fibre space of {£, 11, A^) is 1-dimensional. 
whereas in the rest of the article it is 3-dimensional. Furthermore, ^ is now an element of the 0-sphere {-|-, — } 
and fij = 1 is a constant 0-form. Thus, the analogue to ([^5|) is 

m 

Vv_KI-— iy_F = 0, (31) 

m 



^We use (t, z) rather than (x") to distinguish coordinates on 2- and 4-dimensional spacetimes. 



where V± satisfy the conditions 

g{V+,V+) = -l, g{V+,^^)<0, 

g{V-,V-) = -l, g(F_,^)<0, (32) 

and the only non-trivial Maxwell equation for the 2-form F is 

di^F = -qi.N + qi.N^n (33) 

where iVion = niond/dt is the ion number 2-current and F — Edt A dz where E is the electric field along the 
2;-axis. 

On £, t ^ vl + z^ and the components of the electron number 2-current N = N^d/dt + N^'d/dz corre- 
sponding to pT]) are 



N* = I f{t,z,i,z)dz^a(x+~X^ 



VTTP' 



N^ = I ^==fit,z,t,z)dz = a{^/l + Xl-^/l + Xl], (34) 

where 



[0, Z<X^OTZ>X+ ^ ^ 

with a a positive constant and {X+, X-} scalar fields over spacetime. The 2- velocity fields {V+, y_} satisfy 

and it follows 

N ==a*{vl-Vl). (37) 

Unlike their 4-dimensional analogue, which may include transverse electromagnetic fields, pip are uniquelyo 
solved by 

dV± = —F (38) 

m 

and using 155]) 

_ q2 ^ 

d*dV± = - — (•iV-^Mon) (39) 

m 

subject to the condition dV+ = dV-. 

Alternatively, one may follow the approach adopted in [5] employing a warm fluid model: 



(40) 
(41) 

(42) 



(43) 
where U is the bulk 2-velocity of the electron fluid, 

U^ ^4^^ ^' Z = liV+ + V^), (44) 



ip + p)VuU ^qmuF - 


- Lij{a 


IpA 


g{U,U)^-l, 






.(^,|)<o. 






Here, 






N^nU 







-Proper incorporation of transverse fields requires at least 2 spatial dimensions. 



and, in the electron fluid's rest frame, p is the fluid's energy density and p is the fluid's pressure defined as 
P= \'l + z'^f{t,z,i,z)dz 



P 



VT+P 

It may be shown 



f{t, z, i, i) dz. 



2a 



'1 



n 
2^ 



sinh — 
\2a 



p = ma 



^. 1+f M _sinh-i(^ 
2a\ \2a) \2a 



(45) 
(46) 

(47) 
(48) 



Thus, (|3H [52]) may be replaced by an equivalent field theory expressed in terms of a finite set of moments of 
/ on 2-dimcnsional spacetime. However, the situation is more complicated for waterbags over 4-dimensional 
spacetime where the moment hierarchy is not automatically closed. 

We will now use ([M]) to obtain a non- linear oscillator describing 1-dimensional electrostatic oscillations. 
Let all field components with respect to the laboratory frame (di, dz) be functions of C ~ z ~ vt only (the 
"quasi-static assumption"), where <v < 1, and let [e^^e^) be the basis 



3 = vdz — dt. 



dz — vdt. 



(49) 



The coframe (76^,76^) is an orthonormal basis adapted to observers moving at velocity v along z (i.e observers 
in the "wave frame") where 7 = (1 — u^)^^/^ is the Lorentz factor of such observers relative to the laboratory. 
For example, 7e^(iVion) = — 7'T-ioni' is the ion 1-current in the wave frame. 
In the basis (e-^, e^), V± can be decomposed as 



y± = (M(C) + ^±)ei+V'±(C)e^ 

Note that this is the most general decomposition compatible with equation 
Solving (P2|) for -01 gives 



(50) 
and the quasi-static assumption. 



V'^ = (M + ^±)'-7' (51) 

and additional physical information is needed to fix the sign of ip± . Here, we demand that all electrons described 
by the waterbag are travelling slower than the wave so V'i = ^\/(p + ^±)^ ^ 7^ and ([50|) is 

V± ={^i + A±)e^ - (^i^l + A±f - 7^) '^\^. 

Substituting (|50|1 into equation ([55)) yields 

1 m d/i 
72 g dC' 

and equation p9p yields the nonlinear oscillator equation 

r,2 



(52) 



(53) 



1 d^p 
i^~d^ 



q^ 2 1' 

7 "ion a 

m m 



V(m + ^+)2-7^-v/(m + ^-)^-7' 



with the algebraic constraint 



A, 



A. 



»ion7 •" 

a 



(54) 



(55) 



Longitudinal electrostatic oscillations in 3 spatial dimensions 

We now consider electrostatic waves in 3 spatial dimensions by closely following the above description of 1 
dimensional electric waves. 

To proceed further we seek a form for Wx axisymmetric about i^ whose pointwise dependence in M is on 
the wave's phase C = a;'^ — vx'^ only, where < w < 1. As before, the following results are applicable only if 
the longitudinal component of Vj in the wave frame is negative (no electron described by Wx is moving faster 
along x^ than the wave). 

Decompose V^ in the wave frame as 

Vj = [/i(C) + A{C^)] e^ + V(C\ C) e' + -Rsin(ei) cos{C'')dx^ + Rsm{^^)sm{C'')dx'' (56) 

for < ^^ < TT, < ^^ < 27r where i? > is constant and 



e 



^-vdx^-dx^, e^ = dx^~vdx°. (57) 



Here, (76"'^, 76^, dx^, dx"^) is an orthonormal basis (the wave frame) with 7 = l/Vl — v'^ ■ In the wave frame the 
relativistic energy of Pj = mV^ is m(/x + A)/^ and it follows that y, + A> Q. The component 1/' is determined 
using ([26]) . 



^ = -^[/i + A]2-72[l + i?2sin2(ei)], (58) 

where the negative square root is chosen because no electron is moving faster along x^ than the wave. 
Substituting ([5^ into equation ([77|l leads to 

(59) 



m da r, , 
q d( 




and (HaiiaiSSlEHD yield 




1 d^^l q^ 2 
7^ d(;^ m 


- ^2nR^a 
m 



[P^ + A{e)? - 7'[1 + R^ sin2(ei)]^ sin(ei) cos(ei) d^ (60) 



(c.f. equation ([54]) ) and 

27ri?2 fAie) Me) co,ie)de = -i^^h^ (6i) 

J a 



(c.f. equation ([55]) ) where a is the value of / inside Wx- 

The form of the 2nd order autonomous non-linear ordinary differential equation (|60p for /i is fixed by 
specifying the generator A{e) oi dWx subject to the normalization condition ([6T|) . 

3 Electrostatic wave-breaking 

The form of the integrand in (pD|) ensures that the magnitude of oscillatory solutions to (|60p cannot be arbitrarily 
large. For our model, the wave-breaking value /^wb is the largest fi for which the argument of the square root 
in (pn|) vanishes, 



^wb = maxj - A{e) + -f^l + R^sui^ie) < ^^ < tt j 



(62) 



because /i < ^wb yields an imaginary integrand in ()60p for some ^^. The positive square root in (j62p is chosen 
because, as discussed above, fi + A{e) > and in particular /Xwb + ^(C^) > 0. 

The electric field has only one non-zero component E (in the x^ direction) . Using F = E dx^ A dx^ and (|561 
EZl EH) it follows 



m 1 dfi 

q 72 dC 



E--^^ (63) 



and the wave-breaking limit -Emax is obtained by evaluating the first integral of ()60p between ^wb where E 
vanishes and the equilibriunij value ^eq of fi where i? is at a maximum. Using (|61|) to eliminate a it follows 
that /ieq satisfies 



1/2 



(64) 



with 



since a,v > 0. Equation ([SO]) yields the maximum value i?max of i5, 

V 
Mcq + Mwb+ If X 



(65) 



max ~" ^iiiiiiQYi 



J A{C) sm{£_^') cos{£_^')d£_^' 



Mcq 7T 



[m + Aier - 7'[1 + R^ sin'(e')] sin(ei) cos(Ci)dei d/i 



1/2 



Mwb 



(66) 



Example 



The above may be used to determine a wave-breaking limit for a nearly cold plasma whose distribution's 
transverse extent is much larger than its longitudinal extent. 

Let A{£^^) ~ — acos(^^) where a is a positive constant. Using (IM]) it follows 



E — 9mr)- 
max ^^'^'^lon 



3v 

Mcq + Mwb + ^~ 



Ml-q 1 



1/2 



Mwb -1 

where x = ~ cos(^"'^) and equation ([64|) yields 



(67) 



3w 
2a 



[Mcq + ax] -7'[l + i?'(l-x')] Xrfx=l 



1/2 



(68) 



Equation (|62p may be written 



^wb = max-; - ax + 7a/1 + i?2(l - x^) 
and for a, i?, 7 satisfying 



-1<X<1 , 



a / l + i?2 



i? V a^ + l^R'^ 



< 1, i?> 



(69) 



(70) 



the maximum of /i(x) == ~o-X + 7\/l + i?2(l — x^) over — 1 < x ^ 1 coincides with a turning point x — X 
— cos(^^) of h where 



,,i. a l + R^ 



(71) 



Note that the equihbrium of /i need not coincide with the plasma's thermodynamic equilibrium. 



During the maximuni amplitude oscillation the points S,^ = S,^ catch up with the wave and it follows 



Mwb = -VOTl^W+^fWj. (72) 

H 

For a < i? < 1 equations ([63 EH Ell) yield 

i2^2, ,2 






(73) 



where mcwp •\/2(7 — l)/!?! is the usual relativistic cold plasma wave-breaking limit of E (see, for example, [13]) 
and ujp = y^riionq'^/i'rneo) is the plasma angular frequency. Note that the speed of light c and the permittivity 
So of the vacuum have been restored. The parameter R may be eliminated in favour of an effective transverse 
"temperature" T±cq defined as 

T±cq = 1^ (^cq + Po' ) , (74) 

P:^^ = ma( x^xhx^l (75) 

JWoq 

where Wcq is the support of the distribution with fi ~ fi^q (see footnote [3]) and ks is Boltzmann's constant. It 
follows 



R^J'-h^ (76) 

where the speed of light c has been restored. 

Conclusion 

We have developed a method for investigating the relationship between the shape of a 1-particlc distribution 
and electrostatic non-linear thermal plasma waves near breaking. An approximation to the wave-breaking limit 
of the electric field was obtained for a particular axisymmetric distribution. 
Further analysis of ([66l [64l [62]) will be presented elsewhere. 
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